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Diffusion in multiscale spacetimes 
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We study diffusion processes in anomalous spacetimes regarded as models of quantum geometry. 
Several types of diffusion equation and their solutions are presented and the associated stochastic 
processes are identified. These results are partly based on the literature in probability and percola- 
tion theory but their physical interpretation here is different since they apply to quantum spacetime 
itself. The case of multiscale (in particular, multifractal) spacetimes is then considered through 
a number of examples and the most general spectral-dimension profile of multifractional spaces is 
constructed. 

PACS numbers: 02.50.-r, 05.45.Df, ll.10.Lm, 04.60.-m 



I. INTRODUCTION 

Although many different theories of quantum gravity 
have been proposed to accommodate the gravitational 
interaction and quantum mechanics in a unified frame- 
work, certain features seem to be universal and deeply 
connected with one of the most desired properties all 
these candidates should possess: namely, the absence or 
control of divergences at small scales. One of these fea- 
tures is dimensional reduction, also known as dimensional 
flow QJ-IH . In general, one of the indicators characterizing 
quantum geometry, the spectral dimension ds of space- 
time, changes with the scale, running from ds < 2 (or 
exactly ds = 2) in the ultraviolet (UV) to the usual, 
classical value ds ~ 4 in the infrared (IR). Numerical 
and analytic examples can be found in causal dynami- 
cal triangulations (CDT) [1, [f| , random combs [6l, Ml and 
random multigraphs [1, Q (both sharing some properties 
with CDT), quantum Einstein gravity ( QEG, a lso called 
asymptotic safety^ fiol ITl| , spin foams fl2Tfl5l , Hof ava- 
Lifshitz gravity 
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noncommutative geometry at the 
fundamental [18|,[19| and effective levels |20l - l22l |. field the- 
ory on multifractal spacetimes p3l - l25j (in particular, in 
the realization within multifractional geometry [26H31II) , 
and nonlocal super-renormalizablc quantum gravity |32| — 

H. 

The multifractional framework is rather effective in 
describing geometric and physical features of quantum 
gravity at large. Therefore, it can be regarded either 
as an independent proposal for a fundamental theory or 
an effective framework wherein to better understand the 
multisca le g eometry of the other approaches (examples 
are [2l|, l35l l36j). For this reason, we believe it is im- 
portant to exploit the tools available in multifractional 
spacetimes as much as possible. Due to the young age 
of the proposal, these tools are largely unexplored and 
the purpose of this paper is to continue the investigation 
initiated in Refs. [26H29j . in the meanwhile improving 
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the understanding of dimensional flow as a general phe- 
nomenon of quantum geometry. 

To begin with, we need to refine the discussion on dif- 
fusion processes of Ref . [I?} ■ Systems with anomalous di- 
mension are described by fractional versions of the stan- 
dard diffusion equation, while multiscale systems are ob- 
tained by further but simple extensions. In Ref. [30j . 
which is a condensed exposition of the contents of Sees. 
HVI andlVl the following was pointed out. Once the forms 
of the Laplacian and of the diffusion operator are de- 
termined independently (e.g., by phenomenology), the 
requirement that the solution of the diffusion equation 
be non-negative definite allows one to fix other details 
of the diffusion equation, such as the presence of source 
terms, and to give a probabilistic interpretation to the 
diffusion process. Associating diffusion with a stochastic 
process is an important step towards the understanding 
of the physical properties of the quantum geometry un- 
der scrutiny and of dimensional flow. Multiscale cases 
are more complicated than those with fixed dimension- 
ality, but this does not prevent the construction of sys- 
tems with a composite scale hierarchy. Examples exist 
in the literature of probability and chaos theory that can 
be applied to quantum-gravity scenarios in the present 
interpretation. We shall also construct the profile of the 
spectral dimension for a multifractional spacetime with a 
finite but arbitrary number of characteristic scales. The 
author had the single-scale case in mind when presenting 
the main features of fractal and fractional field theory 
in Refs. [23U261 . [28| . but multiscale extensions also have 
direct applications to quantu m g ravity, like the two-scale 
profile for asymptotic safety (3(J HH . 

The plan of the paper is as follows. Section [XT] in- 
troduces the triple issue of universality, robustness, and 
uniqueness of diffusion and geometric properties in quan- 
tum gravity, and motivates why we expect and study 
nonstandard diffusion equations generated by quantum- 
geometric effects. The main geometric indicators of con- 
tinuum models of fractal spacetimes are reviewed in Sec. 
IIIII with some additional material. After the extension of 
fractional measures to bilateral worlds 12911 . the calcula- 
tion of the spectral dimension of Ref. [271 ] , with further 
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comments on the diffusion probability, can be generalized 
to these cases and to fractional Laplacians. In Sees. MI 51 
MIDI we exploit the recent construction of a momentum 
transform in fractional spaces [29[ and include a class 
of fractional Laplacians not considered in Refs. [U [28j . 
The classification of different types of diffusion is heavily 
based upon known results in the literature on probability 
theory, stochastic processes, and diffusion, although here 
we give the subject a different spin by identifying the dif- 
fusion medium with spacetime itself (Scc. lIVj) . Control of 
multifractional spacetimes is greatly improved in Sec. [Vj 
where dimensional flow is treated analytically. Section 
[VTl is devoted to conclusions. 



II. UNIVERSALITY, ROBUSTNESS, AND 
UNIQUENESS 

To probe the local geometry of a given classical space- 
time, there is a standard recipe founded upon a diffusion 
equation. One places a pointwise test particle on the Eu- 
clidean (i.e., imaginary-time) version of spacetime and 
lets it diffuse from point x to point x' , parametrizing the 
diffusion with an external variable a. It is common to 
assume the diffusion equation }37l - l39| 

(da - V*) P(x, x', a) = , P(x, x', 0) = 8{x - x) , 

(1) 

where the initial condition states the nonextension of the 
probe. The parameter a > acts as an abstract "time" 
variable via the diffusion operator d a , an ordinary first- 
order derivative. The Laplacian V 2 (acting on the vari- 
able x), called more generally the spatial generator in 
probability theory, is a second-order differential opera- 
tor. 

Quantum geometry can emerge either by definition of 
a nonstandard texture of spacetime (as in noncommuta- 
tive and multifractional spaces and in Hofava-Lifshitz 
theory) or from the discretization or quantization of 
the gravitational interaction (as in CDT and asymptotic 
safety), or for both reasons. In all cases, the spectral di- 
mension ds becomes anomalous, i.e., different from the 
topological dimension, and not because of curvature ef- 
fects which exist already at the classical level. Since ds 
stems from the diffusion equation, quantum geometry ef- 
fectively modifies either the diffusion operator d a , or the 
Laplacian V 2 , or the initial condition P(x, x',0), or the 
three of them. To understand the origin of such defor- 
mations of Eq. ([I]) , we also quote some models where this 
happens. 

• Modification of d a . Diffusion "time" acquires an 
anomalous scaling and the diffusion operator be- 
comes fractional (see below), signaling the emer- 
gence of a memory effect along the diffusion flow. 
In the presence of one or more fundamental quan- 
tum scales, the new operator is in fact a sum 
of operators of different orders. To the best of 



our knowledge, there are only two concrete ex- 
amples of a quantum-gravity diffusion equation 
with (multi) fractional diffusion operator: multi- 
fractional spacetimes [l?} (but as an optional con- 
struction) and, perhaps more interestingly for the 
habitues of the field, maybe also quantum Einstein 
gravity (3^|. There, the deformation of da is real- 
ized because the cut-off scale of the theory is not 
identified with the physical momentum as usual 
but, as a powerful alternative, with diffusion time; 
the anomalous scaling is due to the renormalization 
group flow realizing asymptotic safety in the UV, 
while the presence of several scales (and, hence, of 
several diffusion operators in the same equation) is 
guaranteed by the type of action. This liberty in 
the cut-off identification is part of the question of 
universality and robustness outlined below. Also, a 
change in the diffusion operator can be made equiv- 
alent (from the point of view of the asymptotic scal- 
ing of the variance) to a change of spatial genera- 
tor, thanks to a duality between diffusion equations 
which we will introduce later. Finally, the relation 
between number of operators, number of scales, and 
scale-dependence of the spectral dimension is not a 
subject much explored in quantum gravity (but see 
(40| ) and will be extensively analyzed here. 

• Modification of V 2 ,. As the spacetime texture 
gets modified (by quantum-gravity effects or by 
a nonstandard Ansatz of background geometry), 
the differential structure also changes, and with 
that the notion of Laplacian. Noncommutative ge- 
ometry and multifractional spacetimes are two in- 
stances: in both cases, calculus itself and the mea- 
sure of spacetime undergo deep revisions from first 
principles. Qualitatively, asymptotic safety sports 
higher-order Laplacians ( "dual" to lower-order dif- 
fusion operators) because of the anomalous scal- 
ing of the metric within (which survives in inertial 
frames). In Hofava-Lifshitz gravity, higher-order 
Laplacians arise by a simple power-counting argu- 
ment of the quantum theory (l6| . In CDT, spin 
foams, and simplicial gravity in general, the defor- 
mation of the Laplacian is often not explicit but 
one can understand it as coming from the semi- 
classical continuum limit of a discretized geometry 
where multiple scales are introduced [l^]. Dilaton 
black holes are another exam ple of diffusion equa- 
tion with modified Laplacian [411 ] . 

• Modification of P(x,x',0). On a manifold with 
anomalous properties such as in a quantum setting, 
the notion of "point particle" can undergo a change 
to adapt to the different geometric background. In 
other words, the Dirac distribution 8(x — x') no 
longer plays the role of the "delta" and must be re- 
placed by a different distribution. Two cases where 
this happens are quantum manifolds with a min- 
imal length, where pointwise objects are smeared 
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to Gaussians [42J, |43j, and, again, multifractional 
spacetimes, where the Dirac distribution must be 
substituted by its fractional generalization [13, [1^] . 

A classification of the possible diffusion equations should 
help to understand the physics beyond these and other 
models of quantum gravity. To carve concrete inroads 
into multiscale quantum geometry, it is instructive to fo- 
cus one's attention on the characteristics of dimensional 
flow, starting from the most general ones. All models of 
quantum gravity can be equipped with a diffusion equa- 
tion and, hence, a sensible definition of the effective di- 
mension of spacetime (the spectral dimension ds)- Then, 
one can pose a series of questions grouped into three main 
issues: 

Universality: Are there geometric features common to 
all or the great majority of the approaches to quantum 
gravity? If not, are there features which are general to 
many models although not universal? 

Robustness: Are there features of dimensional flow in- 
sensitive to the fine details? 

Uniqueness: Given a model of quantum gravity, is the 
diffusion equation (hence the stochastic process and the 
profile of the spectral dimension) uniquely determined? 
If not, is the physics robust? 

Here we give a qualitative anticipation of what will 
emerge during our analysis. This will help the reader 
in framing the salient characteristics of multiscale and 
multifractal geometries within the mathematical presen- 
tation, placing proper emphasis on those details that play 
a major role in the physics of quantum gravity. 

Part of the introductory motivation was based on the 
observation that a spacetime dimension equal to 2 is uni- 
versally associated with the UV finiteness of quantum 
gravity; dimensional flow from two limiting integer values 
(in particular, ds ~ 2 to ds ~ 4) is an almost universal 
feature of such theories. While the value 2 is itself uni- 
versal, the IR limit ~ 4 is an empirical datum forced into 
the construction of the models, with barely one exception 
26]. 

The number N of asymptotic regimes in dimensional 
flow is not universal because there exist theories, or vari- 
ations within the same theory, predicting more than two 
plateaux in the profile of ds(£), where £ is the probed 
length scale. Asymptotic safety (depending on the op- 
erators assumed in the action) and multifractional mod- 
els (depending on their construction) are two instances. 
CDT show only two regimes but the small-scale one is 
probably an intermediate, not a deep UV feature of the 
model |44| . 

Once N is fixed, how the flow is realized is not universal 
either, inasmuch as different models with the same num- 
ber of characteristic scales generate functionally different 
profiles ds(£). Thus, while the value of each plateau may 
be the same across all possible profiles and the inter- 
polating regimes between them qualitatively similar, the 
actual function realizing the interpolation can be differ- 
ent. Nevertheless, we point out that a certain function 
ds(£) recurs often in the literature and can be regarded 



as general, albeit nonuniversal, for systems with just one 
characteristic scale t*. In four topological dimensions, 
this function is 



d s (£) 



i + (e/i*y 



(2) 



where typically a* = 1/2 and, as we shall comment be- 
low, the power of the length ratio is of little importance. 
The instances where this profile is realized are at least 
the following. 

(a) Multifractional spacetimes. Because this frame- 
work aims to capture, by construction, universal or gen- 
eral multiscale geometries appearing in quantum gravity, 
the derivation of Eq. ()133|) (for isotropic measures) au- 
tomatically contains an explanation of the generality of 
the profile © in single-scale scenarios. It simply stems 
from the most natural adaptation of rods in multiscale 
measurements. 

(b) Dilaton black holes, where Eq. ^ is an accurate 
estimate of an exact but more complicated profile ds(£) 
[4l| . Since the dilaton plays the role of a power- law mea- 
sure weight in the action, technically this case is nothing 
but a radial version of a multifractional spacetime. 

(c) Causal dynamical triangulations. A numerical best 
fit of dimensional flow yields a* w 1/2 and (£/£*) 2 « 
a/50, where a is the dimensionless diffusion time in the 
CDT transport equation Q. 

(d) Random multigraphs, whose geometry was shown 
to be closely akin to (certain approximations of) causal 
dynamical triangulations [H, Q . The analytic form of the 
return probability is, in fact, V{a) = 2G 2 /(2Ga + a 2 ), 
leading to d s = 4[l/2 + cr/(2G)]/[l + a/(2G)], where G 
is Newton's constant. 

(e) Hofava-Lifshitz gravity. There, the scaling of time 
and spatial coordinates is mutually anomalous and gov- 
erned by a parameter z which is fixed in the UV and in 
the IR (to the values 3 and 1, respectively). In general, 
ds = 1 + 3/x (la |: this expression captures the asymp- 
totic behavior of the spectral dimension in the two (UV 
and IR) plateaux corresponding to z = 3 and z = 1. 
The actual profile ds(£) would stem from the exact dif- 
fusion equation featuring (at least two) Laplacians of 
different order. The spatial generator of the diffusion 
equation can thus be written, effectively, as a sum over 
z of Laplacians with z-dependent order [40| . In turn, z 
can be regarded phenomenologically as a scale-dependent 
parameter, which reproduces Eq. ([2]) with a* = 1/2 if 
^) = [l + (^) 2 ]/[(l/3) + (^) 2 ]. 

The latter matching is justified by the evidence we shall 
gather that multiscale (in particular, multi-Laplacian) 
and multifractional diffusion equations yield the same 
spectral dimension profile when the scale hierarchy is mu- 
tually tuned. 

If the number TV — 1 of characteristic or fundamental 
scales of the system is specified, so is the number N > 2 
of asymptotic regimes in the spectral dimension. Profiles 
with different N belong to physically different systems. 
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FIG. 1. Typical single-scale profile of the spectral dimension 
as a function of the probed length scale. (A) Asymptotic 
regimes where ds ~ const (plateaux) and the values of ds 
therein are universal or almost universal. (B) Also interme- 
diate plateaux, possibly reduced to local extrema, are robust 
within a given physical system, but different mathematical re- 
alizations of the same system cannot produce extra plateaux 
or transient features such as bumps, glitches, and so on. (C) 
Details of dimensional flow such as the monotonic slope of the 
profile between different regimes can change with the mathe- 
matical realization but they are physically unimportant. 

Figure Q] shows the typical two- regime (one characteristic 
scale, N = 2) profile. 

The question about uniqueness is more delicate and is 
summarized in the following double claim: 

• Given a physical system, we can associate differ- 
ent diffusion equations realizing very similar di- 
mensional flows with the same physical features, 
i.e., (i) the number N— 1 of plateaux, (ii) their mu- 
tual position, (iii) their values, and (iv) the mono- 
tonicity of the flow between any two plateaux. 

• Other features such as (v) the length of the 
plateaux and (vi) the way they are connected are 
not unique, can vary depending on the mathemat- 
ical implementation of the model, and are not ex- 
pected to play any role as far as physical falsifiabil- 
ity is concerned. 

The number N, the values of ds at the plateaux, and 
their mutual position are almost always default ingre- 
dients, specified by construction either explicitly (e.g., 
requiring that ds ~ 4 in the fR, or introducing only 
N — 1 = 1 fundamental scale such as the Planck scale 
in noncommutative geometries or the simplest multifrac- 
tional setting) or implicitly (e.g., choosing the operators 
in the action in QEG determines N and, in a rather mys- 
terious way, the value and position of the intermediate 
plateau). This should highlight the intrinsically empiri- 
cal nature of the concept of universality. 

Throughout the paper a plethora of examples will show 
that, once N and the values of the plateaux and their mu- 



tual positions are fixed, different diffusion equations will 
be able to generate an analytic dimensional flow connect- 
ing them in about the same way. Obviously, the profiles 
for ds will be different functions of the probed scale, but 
since N is fixed there will be no intermediate local ex- 
trema (maxima or minima, "bumps," "glitches," or other 
transient features like the bump B in Fig. [1} between 
two plateaux. Thus, universal or almost universal fea- 
tures are unique and robust, while the interconnecting 
monotonic slopes in the flow (feature C in Fig. [1} are 
nonunique but robust. 

In transport theory, it is not altogether unusual to asso- 
ciate different diffusion equations (i.e., different transport 
models) with the same physical system and, at the same 
time, being unable to place observational constraints on 
such models. The reason is that, in practice, only anoma- 
lous exponents can be determined by experiments, while 
details of the diffusion process are more elusive. The the- 
ory itself, however, can discriminate among the models, 
because diffusion equations with no probabilistic inter- 
pretation are much less appealing than those with a so- 
lution with stochastic meaning, even if they all realize 
the same anomalous behavior asymptotically. Below we 
will stress and illustrate this point on several occasions. 

In quantum gravity, where we have virtually no ac- 
knowledged observational signature available (with the 
possible exception of the cosmological constant), this lim- 
itation of "uniqueness" is all the more cogent. On top 
of that, there are explanations of why transient slopes 
in the dimensional flow are associated with nonunique 
mathematical features of the models: we refer, in partic- 
ular, to regularization schemes in triangulation settings 
(as in CDT) and in renormalization group flows (as in 
asymptotic safety and Hofava-Lifshitz gravity). There- 
fore, even from the point of view of quantum field theory 
and simplicial geometry it is not always physically mean- 
ingful to expect transient regimes to be unique. Finally, 
we wish to advance the somewhat neglected idea that 
diffusion must have a well-defined stochastic interpreta- 
tion in quantum gravity also. Diffusion, in fact, probes 
geometry itself via a test particle and one should be in 
a position to meaningfully ask what the probability is of 
finding the probe at a certain point. 

III. GEOMETRY OF CLASSICAL 
FRACTIONAL SPACETIMES 

A. Configuration space and HausdorfF dimension 

Continuum models of fractal spacctimcs arc defined via 
an embedding, a measure g, and a Laplace-Beltrami op- 
erator K,. Let M D be Minkowski spacetime in D dimen- 
sions, labeled by Greek indices /i, v, ■ ■ ■ = 0, 1, . . . , D — 1. 
The measure is of a Lebesgue-Stieltjes type [23, [24[ such 
that it can be recast as a Lebesgue measure with a non- 
trivial weight: 

dg(x) = d D xv(x) . (3) 
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We assume that, in momentum units, the scaling dimen- 
sion of the weight v(x) is 



[v(x)]=-D(a-l), 



H A* ' 



(4) 



5 ^ «M < 1 ■ 



V . 



Here the "fractional charges" are parameters asso- 
ciated with each direction. An explicit realization of 
this weight is given by real-order fractional measures 
dg a (x) = d D xv a (x), where (27ll29j| 



*(x) = Y[v a (xn :=n 



(5) 



where the absolute value is taken so that the measure 
is real-valued and positive. This has consequences for 
quantization [3lj . because unilateral measures with sup- 
port only in the first orthant (x M > 0) are eventually 
conjugate to problematic non-negative momenta. Hence, 
we take the integration range to be the whole real line 
for each coordinate. The "isotropic case" = a can be 
assumed, if desired, to illustrate the properties of these 
spaces. Real-order fractional measures make anomalous 
scaling apparent, 



g a {Xx) = \ dn Q a (x) 



Da. 



(6) 



where dn is the Hausdorff dimension of spacetime. dn 
can be computed from a similarity theorem or from its 
operational definition, as the exponent in the power-law 
scaling of the volume of a D-ball with respect to the 
radius, V (£,) (i?) ~ R dH 



B. Laplacians 

The so-called harmonic structure of the fractal 
Minkwoski spacetime M.® is specified by the choice of 
Laplace-Beltrami operator. Consider, for instance, the 
action of a free massless scalar field: 



S 



dg(x) (j)]C(t) . 



(7) 



One can show that a natural Laplace-Beltrami operator 
on these spaces is given by the self-adjoint operator T>: 



1 



, (8) 



where rj^ is the Minkowski metric with signature 
(—,+,••• ,+), Einstein summation convention is used, 
and <9 M = d/dx^ is the ordinary partial derivative. K, 
is a second-order differential operator guaranteeing that, 
upon integrating by parts, the kinetic term in (JT]) can also 
be written in the symmetric form —T>^T>^^> (boundary 
terms vanish on a suitable space of functionals [29j|). 

These properties are highly desirable, or even neces- 
sary, when considering fractal field theory as fundamen- 
tal. As an effective model, however, can be equipped 



with more general non-Hermitian Laplacians. Consider, 
in one dimension, the left and right Caputo derivatives 



UdV)(x) := — — P ^——^/(a/), 

1 (n — 7) J a (x — x'p+ L ™ 



n — 1 < 7 < n , 



(-1)" 



x > a . 



(9) 



dx' 



T{n-i)J x {x'-xp+ 
n — 1 < 7 < n , x < b 



— dim, 



(10) 



where d is the ordinary first-order partial derivative and 
n > 1 is a natural number. A review of the properties 
of Caputo integro-differential operators can be found in 
Ref. [231 and references therein. Here we only recall some 
features. When 7 n, a d n = d n = {-l) n b d n . Left 
and right derivatives are mutually related by a reflection 
x — > a+b—x. In our context, the fractional measure splits 
the real line in two at x = 0, so that the natural fractional 
derivatives are {-^d 1 , +ao B^} (also called the Liouville 
and Weyl derivative, respectively) and {o<9 7 , o<9 7 }, which 
we shall denote, respectively, as 



{oo9 7 ,oo^}, {CM 7 }. 



(11) 



We paired derivatives which are conjugate under reflec- 
tion, in both cases at 0, x — > —x. By conjugate, we mean 



(dy)(x) = (By.)(-x) , 

(9 7 /)(x) = (9V_)(-x), (12a) 

Ud\f)(x) = ( 00 8y-)(-x), 

UBy)(x) = ( oc 9 7 /_)(- a; ), (12b) 

where f—(x) := f(—x). While for ordinary derivatives 
this parity transformation leads to the same operator up 
to a sign, d — > —d, in fractional calculus it transforms 
left into right derivatives, and vice versa. For instance 
@, Eq. (2.3.23)], 



+ 'X 



dx f oa d 1 g 



dxgoad 1 f . 



(13) 



This is the reason why one cannot construct naive self- 
adjoint fractional Laplacians on fractional spaces. In Ref. 
[27| . the fractional operators 



£7 = V^dWl , 



'/ 00 



were defined for a unilateral world where x^ > (a sub- 
script /i denotes the left or right Caputo derivative with 
respect to x^). However, in a bilateral world the coor- 
dinate upon which one differentiates must lie within the 
integration range of the given derivative, so one can use 
both ooi9 7 and oo<9 7 at the same time, or else the pair <9 7 
and d 1 . Therefore, adding some weight factors for later 
convenience, here we shall define the fractional Laplace- 
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Beltrami operatortQ 
1 



^7, ol • — 



y/v{x) 



(14) 



1 



: If 



7 '" Vv(xj 

x (c ^oo9J) k/^) • ,(15) 

where c 7 and c 7 are complex coefficients dependent on 
7, which we choose so that c„ + c n = 1 for n nat- 
ural. Suitable boundary conditions at x = on the 
space of functional yield the properties d 1 d 1 = d 2,1 and 
qiqi = {pi j while for Liouville derivatives it is always 
true that ood^c^ = ^d 2 ^ ood 7 oo^ = oo a r H3, Sec. 
2.3.2]. Thus, the order of the Laplace-Beltrami operators 
([Li]) and (fT5j) is 27, where 7 = £^ 7^/-D. In the isotropic 
case, 7^ = 7 for all /i. In anisotropic formulations, sum- 
mation with the Minkowski metric must be accompanied 
by suitable dimensionful coefficients, which we omitted 
in Eqs. flU]) and (JTSJ. When 7=1, £i )Q = /C = K 1>a ; 
when 7 = n and a = 1, /C n ,i = K, n ,i = rj^d^d™. In 
Euclidean signature, we denote the ordinary Laplacian 



In a field-theory setting, combinations of Laplacians 
with integer 7 > 1 will, in general, introduce ghost modes 
(e.g., [46M48j ) . while fractional Laplacians (0 < 7 ^ N) 
are typically associated with continuum quasiparticle 
spectra. The first fact already foretells an issue with 
values 7 > 1 which will be further strengthened when 
looking at probabilities in diffusion processes. 

For later use, we recall that (e.g., [27| ) 



<9 7 £ 7 (Ax 7 ) = A£ 7 (Ax 7 ), 

3V a = Xx 1 ' 1 Ei, 2 -y{Xx) . 

where E is the Mittag-Lcffler function 

E a {z) :=E aA (z), 

E a ,b( Z ) : = ~ ~' a>0 



T(an + b) ' 



(20) 
(21) 



(22a) 
(22b) 



Analogous expressions hold for the right derivative with 
x — > —x. Also, 



,d~<e Ax = (-A) 7 e 



(23) 



^ d 2 
72 _ ° 



(16) 



and, in general, higher-order differential operators as 



T — 



n G N. 



(17) 



With the measure ((5]) and bilateral integration, from 
Eq. (|13|) one sees that the operator (fl~5|) (with coefficients 
chosen as below) is self-adjoint with respect to the scalar 
product of the Hilbert space of suitably "good" functions 
given by the intersection of the domain of /C 7 Q , and the 
space of functions [29[ which are L 2 with respect to the 
measure g. Given two such nontrivial functions / and g, 
the scalar product is defined as 



if, 9) 



dg(x) f*(x)g(x), 



(18) 



where the asterisk indicates complex conjugation, so that 

(/ , £7,00) = (£ 7 ,a/ ' 9) & C 7 = C*. (19) 



c* in the definition ([Toll . This 



From now on, we set c 7 — 
argument is sufficient to exclude /C 7jQ! (which cannot be 
Hermitian on this space) but we shall give another moti- 
vation soon. 



1 We call them "fractional" even if in many formulas the measure 
v(x) may be more general than a fractional measure v a (x). 



C. Momentum transform 

We restrict now the measure to be non-negative def- 
inite and factorizable with respect to its coordinate de- 
pendence: 

v(x) = l[v fi (x>*), «„(a^)>0, V^, (24) 



where the index /1 in is not vectorial. The fractional 
measure ([5]) satisfies these conditions. Let Ai® be an- 
other space in a Z3-dimensional embedding, spanned by 
"momenta" k and with measure 



dr(k) = d D kw{k) 



Da' 



(25) 



where a' is not necessarily equal to a. If w(k) obeys 
the same properties as v [as in the case of two fractional 
measures v = v a (x) and w = v a i(k)], then there exists 
a family of unitary invertible mappings F : M.® — > 
identifying as momentum space [29(. This is the gen- 
eralization of the Fourier transform on Lebesgue-Sticlt jcs 
spacetimes with factorizable measures: 



+00 



dg(x)f(x)<B*(k,x) =: F fllT [/(x)],(26a) 



f(x) = / dT(k)f(k)<B(k, X ) 



(26b) 



where <e(fc, x) are the kernel functions on which the trans- 
form is expanded. By construction, they give an integral 
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representation of the fractal Dirac distributions in posi- 
tion and momentum space: 



dr(k) <b* (k, x)s(k, x') — 5 v (x,x) 

S(x-x') 



\J v{x)v(x') 



+00 



dg(x) <s*(k, x)<s(k' , x) — S w (k,k') 

S(k-k') 



w(k)w(k') 



, (27) 



, (28) 



which are quite generally nontranslationally invariant. 
When v — w, F g ^ T is an automorphism, the 5 distri- 
bution is the same in both spaces and the kernels e are 
symmetric in x and k. 



D. Eigenvalue equations 

In Ref. [29 ]. it was shown that the <b's are either 
weighted Bessel functions of the first kind or weighted 
phases, where the weight is [w(k)v(x)]^ 1 ^ 2 . For calcu- 
lational purposes, we demand that the kernel functions 
<b are eigenfunctions of the chosen Laplacian operators. 
However, from Eq. (|20[) one sees that none of the kernels 
can be eigenfunctions of the operator /S 7jQ , so we cannot 
use the momentum transforms to diagonalize this Lapla- 
cian. From now on we drop JC 1>a from the discussion. On 
the other hand, phases are eigenfunctions of the ordinary 
derivative, of the left derivative from —00 to x and of the 
right derivative from x to +00 [Eq. ([23]) ] . In particular, 
the phaselike kernel 



e(fc, ; 



y/w(k)v(x) (2tt)* 



(29) 



is an eigenfunction of both K, and /C 7jQ! : 



K.<R(k,x) = —k 2 (B(k,x), 

k '. k^k^ — — - /lq ~\~ k^ ~\~ ■ 
/C 7:Q , <e(fc, cc) = F 2j (k) ce(fc, x) . 



(30) 



k 2 



(31) 



We now determine the form of the eigenvalue F 2j (k) = 
F 2 ^{k [) ,k 1 ,...,k D ~ 1 ). Take a function f(x) = h(kx), 
where k 7^ is a constant (later one can analytically 
continue the final result to k = 0). From Eq. ([TU|) one 



finds that 

G»#7)(a0 = {oodf h)(kx) 

z—kx 1 f . 1 

= — / dz 

T(n - 27) y_ sgn ( fc )cx> 

k n ~ Y 

X [fc-i(z-z')] 27+1 "" 9 "' /l(Z ' ) 

= k 2 ^e(k)Udph)(z) 

+ (-k) 2 -<6{-k) (ood?h)(z) 

= |fc| 2 W)(oo^)(z) 

+e(-k)( 00 Bph)(z)}, (32) 

where 8 is the Heaviside distribution: 



0(k) 



1, k > 0, 
0, k < 0. 



(33) 



In what follows, the value of 6 at k = is irrelevant as 
7 > 0. Similarly, for the right Liouville derivative 

U^/Xs) = \k\^[6(k) (ooB?h)( Z ) 

+6(-k)( 00 d^h)(z)], (34) 

so that we can combine Eqs. (|3"2"j) and ([M)) into the op- 
erator 

(A 7 f)(x) := c 7 Ud^f^x) + c 7 Ud 2 x V){x) 

= \k\ 2 ^(k)Udph)(z) 

+c 1 {-k){ 00 d^h){z)], (35) 

where c 7 (fc) := 9(k) c 7 + 0{—k) c 7 . 

For the function f(x) = h(kx) = e Afc:c , we get 



A- 



[c 7 (fc) A 27 + c 7 (-fc) (-A)^]|fc| 2 V fa , (36) 



so that, for A = ±i = e ±l '*'/ 2 , we have 



A 7 e 



[c 7 (fc) 

"sgn(fe)7 



+ c 7 (-fc)e T " 7 pl 



\k\ 2l e ±ik3 



(37) 



Some remarks about the eigenvalues of A 7 follow: 



(a) Choosing c 7 7^ c* can lead to complex-valued spec- 



tra. 



(b) Otherwise, for c 7 = c* the spectrum of A 7 has 
a discontinuity at k = and, thanks to the self- 
adjointness of the Laplacian, it is real. Writing 



c 7 = |c 7 | e l7Tip , one has 



"sgn(fc)7 



ct := 2|c 7 |cos[7r((^±7)] 



(38) 



(c) If c 7 



l/(2i) = -1/2), then c± = 



± sgn(fc) sin(7T7) and the signature of the spectrum 
depends on the sign of both k and the phase. 

(d) If c 7 = 1/2 (<p = 0), then = cos(-7T7) depends 
neither on the sign of k nor on that of the phase, 
and the spectrum is semibounded. When 7 = n + 
1/2 is a half integer, n £ Z, plane waves are in the 
kernel of A 7 . When 7 = n, c 7 = (—1)™ and the 
spectrum is positive (negative) definite for n even 
(odd), with zero eigenvalue at k = 0. 
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Moving to many dimensions and the operator /C 7jCt , we 
conclude that the eigenvalue in Eq. (|3Tj) is 

F 2 "(k) ■= -4 en(k o h |fc°| 27 + £ <n (fc . )7 • (39) 

i=l 

Since the sign of the phase should not play any role, we 
will often choose the coefficients 

1 

c 7 = - e l7r(1 - 7) , F 2 ~<(k) := |fc°| 27 - l^l 27 ' ( 4 °) 

i=l 

Other operators are possible whose eigenvalues have 
fractional momentum dimension. Rotation-invariant ex- 
amples are Riesz Laplacian I45l . Sees. 2.10 and 5.5] and 
the Riesz-Bessel operator [49l [50j. and, in Lorentzian 
signature, the Lorentz-covariant fractional powers of the 
d'Alembertian, (□ — to 2 ) 7 , corresponding to a spectrum 
F 2 ^{k) = -(fc 2 + m 2 )T 

IV. DIFFUSION EQUATIONS AND 
STOCHASTIC PROCESSES 

The s pec tral dimension of M.® is defined via a diffusion 
process [23] . Here we extend that calculation to the case 
of bilateral measures and fractional Laplacians, as well as 
to various classes of stochastic processes associated with 
different ranges in the diffusion exponents (3 and 7. We 
review these processes from the literature of stochastic 
theory and comment on the physical meaning of the heat 
kernel as diffusion probability. 

To simplify the logic presentation, it is convenient to 
specialize first to the case of fixed dimensionality (no 
scale hierarchy). We concentrate on the continuum for- 
mulation of fractional calculus, which guarantees anoma- 
lous (in particular, fractal) geometric properties of space- 
times [27H29f and anomalous correlations in diffusion 
problems (e.g., (5814621 ]). Another assumption we make 
here is to ignore curvature. The latter modifies the spec- 
tral properties of spacetime even in a classical setting, 
except in the UV limit a — >• 0. In the absence of di- 
mensional flow, one can include this limit in the defini- 
tion of the spectral dimension, but in general this would 
prevent the probing of large scales. Effective quantum 
spacetimes, anyway, are modified by quantum geometric 
effects even in the absence of curvature, both globally 
and locally in inertial frames, which motivates the as- 
sumption^ 



2 Another point of view is that a multiscalc geometry is very of- 
ten probed in "snapshots" taken at different given scales; these 
snapshots are then collected to give a fragmentary picture of di- 
mensional flow, typically in the UV and IR limits and in transient 
regimes. Each snapshot is related to a fixed spectral dimension 
and the limit a — > is legitimate. In our approach we see global 
modifications of spacetime and we have full control of dimen- 
sional flow, so we do not adopt this perspective. 



As in the ordinary case (JIJ , the geometry is probed by 
a process governed by a hyperbolic diffusion equation (or 
fractional wave equation, for 1 < /3 < 2), but of the form 

(d? - K* >a )P(x, x',a) = S(x, x',a), (41a) 
P(x,x',0) = 5 v (x,x'), (41b) 

where d@ is the Caputo diffusion operator (we do not 
consider 00 3^ since is the natural derivative for a pro- 
cess with memory loss starting at a = 0), K^ a is the 
Euclideanized version of Eq. ([T5|) . P(x,x',cr) is the heat 
kernel and S is a source term. The external "time" or 
"scale" a has actually dimension [a] = —2~f/f3. As far as 
the author knows, the source S has been set to zero in 
all diffusion processes in quantum gravity considered in 
the literature. 

Different ranges in the parameters /? and 7 can corre- 
spond to physically inequivalent diffusion processes, for 
which there exists extensive literature (for a = 1 and 
other types of fractional Laplacians) [58Tj6ll . [63| . Set, 
first, a = 1 and 5 = 0. 

(a) When /3 = 7 = 1 (Eq. ((TJ), one has ordinary diffu- 
sion with the solution 

2 

P(x,x',a) = in(r,<r) := 6 * , (42) 

(47T(j) 2 

where r 2 := ^ \x^ — x'^\ 2 is the Euclidean dis- 
tance between the two points. The diffusion equa- 
tion is associated with a Wiener process B(o~), also 
known as standard Brownian motion. A Wiener 
process is such that (i) B is continuous in a al- 
most surely (i.e., with probability 1), (ii) 5(0) = x' , 
and (iii) the increments of B are independent and 
governed by the Gaussian distribution m, so that 
B(o-)-B(a') ~ui(0,o--cr / ) for a' < a. Here B is 
the random variable denoting the position x of the 
particle at time a (in D dimensions, it is a vector). 

(b) When < /3 < 1 and 7=1, the process possesses 
a heavy tail in waiting times, leading to a delay of 
particle diffusion and, hence, to subdiffusion^ 

(c) When j3 = 1 and < 7 < 1, the process is a 
Levy process. A Levy process L has the follow- 
ing properties: (i) L(a) is right continuous with 
left limits almost surely, (ii) L(0) = x' almost 
surely, (iii) for any sequence a n -i < o~ n , the in- 
crements L(o~ n ) — L(a n —i) are independent, and 
(iv) their distribution is equal to the distribution 
of L(a n — cr„_i). It is a self-similar process like 



3 The fractional-time diffusion equation is not associated with frac- 
tional Brownian motion \&i |65ll , whose probability density 
function is different from the solutions in Sec. IIV Dl 
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fractional Brownian motion, but with an island- 
type structure, characterized by a heavy-tailed dis- 
tribution and "long jumps" connecting clusters of 
shorter steps. This heavy tail in space is associated 
with superdiffusion. 

(d) Systems outside the range 

0</3,7<l (43) 

correspond to transport pseudoprocesses for which 
Eq. (|4ip may not provide a well-defined diffusion 
model in a probabilistic sense (i.e., the signature 
of P may be indefinite, depending on the case and 
on the presence of a source <S ^ 0). Probabilistic 
processes do exist for 1 < f3 < 2 when 7 = 1 or 
7 > /3/2 in one dimension [63]. Models of other 
systems are still unknown, but there exists a wealth 
of studies on pseudoprocesses. 

For a 7^ f, one introduces extra friction and potential 
terms. When (3 = 7 = 1, diffusion is said to be normal 
or Gaussian (the mean squared displacement of the test 
particle grows linearly with diffusion time), otherwise it 
is anomalous. 

In general, processes described by a fractional or 
higher-order diffusion equation arc not Markovian, mean- 
ing that, in contrast to Markov processes, future states 
depend on the present state but also on past ones, and 
the Chapman-Kolmogorov equation docs not hold. This 
feature is expected in the case of fractional diffusion equa- 
tions, where memory drag is typical of pseudodifferential 
operators. However, it holds also for integer higher-order 
equations where the initial condition appears explicitly 
as a source term. 



A. Solution scheme and probabilistic 
interpretation for S — 



where 

F 2l {k) = -^Ifc^l 27 (48) 

is the Euclidean eigenvalue of the Laplacian with the 
"Riesz choice." Notice that, for 7 = a, —F(k) = \k\ is 
the 2a- norm of fc M , naturally equipping fractional spaces 
[13. Explicity, Eq. (JUJ) reads 

P(x, x',a) = 1 =up(x - x' ,a) , (49) 

\J v(x)v(x') 

r+00 ,Du 

u p {x-x\a):= J /fc (a)e ifc -(*-*). (50) 

For later convenience, we make the [3 dependence explicit 
via a subscript in u. The measure prefactor shows that 
P is not translationally invariant. 

The full calculation of P is not actually needed to 
get the spectral dimension, but it is extremely useful 
for clarifying the physical meaning of the diffusion equa- 
tion. When the Riesz Laplacian is chosen, the integral 
can be done exactly, due to the fact that fk is rota- 
tionally invariant. However, for the fractional Laplacian 
(|T5|) with eigenvalues (|48|) the solution of the multidimen- 
sional fractional diffusion equation is considerably more 
demanding than in the one-dimensional case@ Here, we 
give an incomplete discussion of the properties of P. 
First, we argue that, in D dimensions, P is not positive 
definite for 7 > 1. Then we find the exact solution and 
its asymptotic behavior in D = 1, and then in D generic 
for (3 = 1 or 7 = 1 . Many of the steps are adaptations 
of their analogs in ordinary spacetimes to the case of a 
fractional textureH 

For general (3 and 7, fk is the Mittag-Leffier function 



To solve the diffusion equation, assume v(x) is factor- 
izablc and write P and S in terms of the kernel functions, 

/+00 
dT{k)f k (0-)<B(k,x)<B*(k,x'), (44) 
-00 
r+00 

£>(a;, a/, <t) = / dr{k) Sfc(<r) <s(k, x)<B*(k 7 x 1 ) , (45) 



where the function fk must solve the fractional differen- 
tial equation 

[dP - F 21 (k)]f k (a) = s k (a) , f k (0) = 1. (46) 

Notice that the initial condition in Eq. (|4"Tj) is respected, 
by virtue of Eq. ([27]). 

Setting Sfc = for simplicity, the solution of Eq. (|46j) 

is 



fk(o-) 



TF 2 -<(k) 



/3 = I 



Ep[F 2 ~<(k)o f3 ], (3 general, 



(47) 



4 Its form is akin, but not equal, to the heat kernel in D di- 
mensions, for which several formulations and techniques (ex- 
act, asymptotic, and numerical) have been developed. For sim- 
plicity we consider the case where the fractional order 7 and 
the fractional charge a are the same in all directions, although 
anisotropic configurations 7 M ^ 7„ can be of much interest [66l - 
[7lll . Anisotropic transport is, in fact, typical of Hamiltonian sys- 
tems with noncrgodic dynamics and a hierarchical set of islands, 
where particle jump distributions are asymmetric. Examples in 
nature are porous media in geologic aquifers \TA [73ll . 

5 A possible source of confusion might be that the fractional diffu- 
sion equation with a = 1, i.e., the one considered in the literature 
of chaos theory and percolation systems, already describes a dif- 
fusive process on a fractal structure. Analogously, in quantum 
gravity one could interpret the fractional wave equation with 
a = 1 as a diffusion equation on a fractal spacctime; this is 
the case of QEG, for instance, where a = 1 and spacctime does 
have fractal properties. The introduction of a nontrivial measure 
weight v(x) in position space further changes the background 
geometry and topology, as happened in some first attempts to 
generalize the diffusion equation to fractals [7347(1 . 
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of order /3. Equation (|49j) implies the normalization 



d D xv(x) J — — —P(x, x' , cr) 



v(x) 



+00 



d D x \J v{x)v{x') P(x, x , a) 



—00 
+00 



d xup{x — x' , cr) = 1 , 



(51a) 
(51b) 



where we chose the natural symmetrization in x and 
x' . The extra factor in Eq. (|51l) suggests to con- 
sider P(x,x',a) = ^/v(x')/v(x)P(x, x' , a) = up(x — 
x',a)/v(x) as the probability density defining the frac- 
tional system. The extra factor yj v(x') /v(x) would leave 
both the initial condition and the return probability un- 
changed (it is 1, respectively, on the support of the delta 
and when x = x') but it would imply that the Laplacian 
in Eq. (|41|) does not correspond to the self-adjoint oper- 
ator appearing in the dynamics. In particular, it should 
be of the form JC = v~ 1 (x)d fl d' 1 [v(x) ■ }. This can be jus- 
tified by unravelling the natural structure of the proba- 
bility density as a fractional bilinear; this point will be 
further discussed elsewhere. 

If we interpret the heat kernel as the probability to find 
the particle at x at time a with initial position x' at time 
(7 = 0, then we must require that up > 0. Consequently, 
we can define the moments along the direction fi as 
+00 



{\x-x'\ a ) ll := / d D xup{x-x\a)\x ll -x' IJ \ a , 



a > 0, 



(52) 



where we expressed the integral already in terms of up 
since measure factors cancel one another. The defini- 
tion of the moments allows us to show that up (hence 
P and P) is not positive definite for 7 > 1. In gen- 
eral, not all moments will be finite but we assume that 
(\x — x'\ 2 ) ^ < 00. This is the case in D = 1 and for all the 
D-dimcnsional special cases below; it is quite reasonable 
to expect finiteness of the second moment in general, lest 
the walk dimension be ill defined. Assuming also that 



up > 0, we have (\x 
momentum space, 



> 0. Transforming up to 



j-D / 1 \ ik-(x—x ) 

a xup(x — x ,a) e y 1 



(53) 



one also notices that 



=0 



k u =0 



yy nn+1 27w(27w-l)^ 
^ ' T(0n + 1) 



fc u =0 



This expression identically vanishes if 7 > 1, which con- 
tradicts the assumptions. Therefore, P is not positive 



definite if 7 > 1. This conclusion holds only if S = and 
it cannot say anything about cases where the diffusion 
equation has a source. 



B. General solution in one dimension for 5 = 

In D = 1 dimension, a = 1 and in the absence of 
source, one can find the exact solution of 



(d^-d^)up(x-x',a) = 



(54) 



and prove that up is non- negative definite for the range 
(|4"3"|) . The solution is the same as for the Riesz Laplacian 
(e.g., [Hsl, [HI), and in fact we denoted as dj 2 ^ the Riesz- 
type spatial generator, i.e., the operator /C 7j i in one di- 
mension. The integral to solve is the Fourier transform 
of the Mittag-Lefflcr function. Rescaling the coordinates 
to the dimensionless variables 



J. , 



a 2 ~> x , 



(55) 



we get 



+00 



dk 
2^ 



Ep{~\k\ 2 ^a p )e lk( - x '- x) 



(-1)" 



n=0 Vr ; 



+00 



dk 



2n 



2771 ik(x — x') 



+00 



r(2 7 n+l) 



71=0 

X cos 



_£/|2 7 n+l r(/3n + 1) 
[|(2 7 n+l)] , 



where we used formula 17.21.25 of Rcf. [77|, valid for 
nonintcger 7. The series is absolutely convergent if 27 < 
P and can be resummed exactly for B = 27 and also for 
certain values outside this range |7£|. The (translation- 
invariant part of the) heat kernel is thus 



up(x — x 1 ,<t) 



+00 

1 



, -,y 1 +i sin ( 7r 7^) r(2 7 n+ 1) 

1 ' 7T r(^n + i) 



x/|2 7 n+l 



(56) 



Positiveness of up (and, hence, of P) is guaranteed if the 
fractional exponents are in the range (|43p (in our case 
we have also a measure prefactor but it is positive); the 
proof can be found in Ref. [78| . 

For \x — x'\ 3> 1, the asymptotic limit of up is given 
by the n = 1 term, 



up(x — x' , cr) 



\x-x'\ 2 ^>a /3 ^ 



sin(7T7) T(27 + 1) 



r(/J + i) 



IX — X 



'12-7 + 1 



(57) 
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FIG. 2. The probability densities (|58H (dashed curve) and 
P = up/v a (x) (solid curve) with f3 = 1/2 = a, x' = 1 and 
(j=l. The plot of P = Uf}/ ^/va (x)v a (x') is very similar to 
the one for P. 



This expression implies that all moments (\x — x'\ a ) are 
finite if < a < 2-f < 2. 

When 7=1, the above expressions are ill defined and 
one should find the solution in a different way. Using the 
Laplace-transform method, one can show (451, Corollary 
6.5] that, if < j3 < 1, the solution is 



up(x — x', a) 
where 



a 2 



W[-^l-^-\x-x'\) , (58) 



+ 00 



W(a;b;z) := £ 



n=0 



n\T(an + b) 



(59) 



is Wright's function. An exact solution also exists for 
1 < /3 < 2, which coincides with the one above if 
d a Ufj(x — x',0-)|<j = o = is imposed [III Corollary 6.6]. 
The solution (|58[) and its fractional-space counterpart are 
plotted in Fig. [2J 



C. Levy process (<3 = 1, < 7 < 1, 5 = 0) 

The diffusion equation associated with Levy processes 
is integer in diffusion time and fractional in the spatial 
generator (i.e., the Laplacian). Our Laplacian is not the 
Riesz operator as in standard Levy processes (except in 
D = 1) but the physics is qualitatively the same. The 
diffusion equation is then 




/ 1 Nn+i sin ( 7r 7^)r(27n+l) 
[ ' 7T " r(n + l) 



(61) 



To get an asymptotic expression for x ~ x' , we notice 
[791 ] that the integrand in the momentum expression for 
u\ is even and one can take the real part of e lki - x ~ x \ 
Expanding the cosine, we obtain in one dimension 



+00 



e -cr\k\ 2y e ik(x'-x) 



dk 
2^ 



J2^\ } (x-x') 
' 7r(2m) 

m=0 v ; 



+00 



(-1)' 



^ (2m)! 2tt7 



2m+l 

^ > a -^ {x -x<r\ 



where we used [77J, Eq. 3.478.1]. Then, 

D 

Ui(x,x',a) = Y\_ ui(x^ - x'^a) 

Vl) ? 




r 



( 2m+l \ 
\ ^ ) 



(2m)! 2tt7 



a y (x M - x' ) 




(62) 



Equation can be modified to the case where the spa- 
tial generator is the Riesz-Bessel operator (49[. 



D. Fractional-time diffusion equation (0 < /? < 1, 

7=1, 5 = 0) 

Some subdiffusivc processes are characterized by a 
transport equation with fractional diffusion operator and 
a second-order ordinary Laplacian [53, [63, l80H84| 



V*) u =O. 



(63) 



A nontrivial measure dependence will be present in the 
fractional-space case, but this generalization will be im- 
mediate. When 7=1, the integral in Eq. (|49|) can be 
factorized into a radial and an angular part. Eventually, 



(8 - K* a ) ui = . (60) 

When p = 1, Eq. (gTJ) is the product of D Levy dis- 
tributions <s I fe I ~" with Levy index 27 [59|, Sec. 4]. The 
expression of the heat kernel is (D copies of) Eq. (|56l) 



6 Often it is also described by another transport equation, called 
a bi-fractional or fractional Fick equation, where the diffusion 
equation is "redistributed," (d a - dl~ d%)P = [H, IH. For 
Caputo derivatives these two formulations coincide, while for the 
Riemann-Liouville derivative a source term must be added to the 
first equation. 



12 



up can be cast as a Fox function H 2 ' 2 for < (5 < 1 [451 . 
Corollary 6.7], or as a Meijer G-function as we shall do 
later; if one assumes that d a up(x — x',cr)L=o = 0, the 
same solution holds also for 1 < j3 < 2 45, Corollary 
6.8]. up (hence, P) is nonnegative definite for all these 
values of /?, including /} = 2 [dS. [85|. 

Another formulation, more convenient to study the an- 
alytic properties of P, is the following [86] . In one dimen- 
sion, if < j3 < 1, the solution of 

(dP-d%)up(x-x',a)=0, (64a) 
up(x-x',0)=S(x-x'), 0</3<l, (64b) 




up(x — x', a) = 



+ OG 



ds 



e 4 " 



u 2 p(x-x',a), (65) 



where s is a parameter with engineering dimension [s] = 
[er]/2 = — 1//3 and u 2 /3 is the solution of the analogous 
problem 



{d 2 f -dl)u 20 (x-x',a) = O, 
with initial condition 

u 2 p(x-x',0) = 6{x-x'), 0<f5<\ 

or 

f*2/3(iC — 0) = S(x — x') , 



(66) 
(67) 



d a u 2 p{x-x',a)\ a=0 =0, i < /3 < 1 . (68b) 



Clearly, the initial condition in Eq. (|64|) is respected. In 
particular, when (3 = 1/2 the solution is the integral 
product of two Gaussians, 



Mi (x — x' , cr) 



-t-oo 



(x-g') 2 



f/.S 



e 4" e 



(69) 



and the generalization to D dimensions is straightfor- 
ward. The distribution 



ui(r, cr) 



e 4<r e 4s 



V 7TCT (47TS)' 



(70) 



is the solution of the diffusion equation 

(dl - \?l)ui(x-x',a) =0, (71a) 

!ii(i-s',0) = 5(x-a;'), < /? < 1 . (71b) 

All the above solutions are manifestly non-negative defi- 
nite and normalized, 



up > , 



d x Ufi(x — x\ cr) = 1 



(72) 



so they can be interpreted as probability densities. Under 



-ft/2 



X, s 



the rescaling to dimensionlcss variables x = a 
a~^l 2 s, the cr dependence of up can be factorized as 

up oc cr 2 . (73) 



FIG. 3. The solution (|74|l of the quartic diffusion equation 
(|84[) (/3 = 1/2, dashed curve) and P = u 1 / 2 /v a (x) (solid 
curve) with D = 4, a — 1/2, a/ = 1 and a = 1. P and U\/2 
are non-negative definite and can be interpreted as probability 
distributions. 



The integral in Eq. (|70j) can be done exactly to yield a 
Meijer G function: 



ui(r, cr) 



1 



Lr, 



:S() 



/ r 4 










\256<t 





1 

2 


'0 



(74) 



which is positive definite (Fig. [3j). 

When 7=1, Eq. f[57j) vanishes and all moments with 
a > arc finite. An explicit determination of the mo- 
ments for even a [78j makes use of the diffusion equation. 
In fact, 



-oo 



-oo 



/ \a 



^(Ix-x'Hm = / d D x\x^-x' fl \ a d^(x-x',a) 

) 

d D x\x fl — a;JJ a K.® s u i (x — x' , cr) 
d D xui(x — x',a) JC^^x^ — .t^| 

J — oo 

so that for 7 = 1 

/+oo 
d D xu i (x-x',*)d?\x fi -x^\ a 
-oo 

/•+oo 

= a(a — 1) / d D xui{x — x' ,a) 

(75) 



Using the fact that (e.g., H3) 

we obtain the second moment (a = 2) 

2 



(|s-*T) 



M r(/3 + i) 



(76) 



(77) 
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and, by recursion, all even moments: 

(\*-A 2n )» = l? a n + l \ ° \ neN. (78) 



Odd moments vanish, (\x — x'\ 2n+1 )^ = 0. In particular, 
the average of the square Euclidean distance r 2 covered 
by the particle is a function of diffusion timeO] 



<(r-(r)) 2 ) = 



2D 



r(/3 + i) 



(79) 



E. Iterated Brownian motion (/3 = 1, 7 = 2, <S 7^ 0) 

An example of particular i mportan ce for us is iter- 
ated Brownian motion (IBM) [8a4l04 |. To illustrate it- 
erated Brownian motion, we set a = 1 and consider 
the ordinary Z?-dimensional higher-order operator V n , 
Eq. (|17p. Given a process X(a), let E be the expec- 
tation associated with X(0) = x' . Then, the function 
u(x—x\ a) = E{f[X{a)]} solves a certain diffusion equa- 
tion with initial condition u(x — x' , 0), where x' is fixed. 
For an ordinary Brownian motion X(o~) = B(a), the dis- 
tribution u obeys the diffusion equation 

(d a - V 2 ) u(x-x',a) =0, u(x-x',0) = S(x-x') . 

(80) 

Let B± be two independent standard Brownian motions 
defining the two-sided Brownian process 



X(t) 



B+{t), t>0, 
B_(-t), i<0, 



(81) 



where ( e I. Let B be another independent Brownian 
motion. The IBM is the process defined as 



X mu (cr):=X[B(a)}, a > . 



(82) 



where the Brownian motion B acts as a clock to the two- 
sided motion. For this reason, iterated Brownian motion 
is often called also Brownian-timc Brownian motion. An- 
other definition of IBM is, given two independent Brow- 
nian processes B12, 

X{ BM (a):=B 1 [\B 2 (a)\], a > . (83) 

One can extend these definitions to multiple iterations of 
n processes. If X = I™ 01 X = Aj' BM , the diffusion 
equation is ^BM, HMUqI 



(<9<j — V x ) u(x — x , a) 



V 2 x u{x-x',0) (84) 



7 In fractional spaces, the natural distance is a 2a-norm |27| but 
here, for simplicity, we use a (2ap)-norm with p = 1/a > 1, 
which is topologically equivalent. Restoring measure factors and 
extending the discussion to fractional spaces does not entail fur- 
ther difficulties. 



in D dimensions. Since the process is non-Markovian, the 
initial condition u(x — x',0) explicitly acts as a source. 
Other processes than Brownian motions can be composed 
toge ther and generate different diffusion equations (see 



t 104| for an account). In Se c. IVBIwe shall consider tele- 
graph processes [H [H, [§1 1105L Il06l | . 



A remarkable duality shows that fractional and it- 
erated Brownian motions can be identified. Namely, 
Eq. ((701) is a solution of both ((711) and f84]). In gen- 
eral, there exists a triple connection between fractional 
diffusion equations with fractional time <7, higher-order 
diffusion equatio ns with integer time, and iterated pro- 
cesses [H M, M, DM Ho2l . fio4| . 

When the kinetic term is a second-order differential 
operator, the generator of the IBM is a "half derivative" 
which can be rigorously defined [95j . Intuitively, this 
stems from the iteration of the fractional diffusion equa- 
tion (|7"T|) . Let U1/2 be the solution of Eq. ([TT]) . Applying 

da twice and using the property d 1 / 2 d 1 / 2 = d, valid for 
Cap uto derivatives [see, e.g., Eqs. (2.50)-(2.53) of Ref. 
[23], we have for x - x' ^ 



cLit 1 



= 83 dS ui 
^ dhlu 



(85) 



For x ~ x' , the third step is not correct and one can 
check that 1/4/2 is indeed a solution of ([84")) in a weak 
sense (i.e., by integrating the diffusion equation with a 
test function); this is the point where the source arises. 
The interested reader can find the proof in the references 
cited above. 

The equivalence between fractional and iterated Brow- 
nian motion can play a role in the interpretation of quan- 
tum geometry at the UV fixed point, including QEG j36| . 
For this reason, it is important to stress the physical 
meaning of transport equations such as (|7Tj) and 
There do exist physical systems which are associated with 
these hyperbolic equations. In fact, iterated Brownian 
motion provides one stochas tic d escription (among oth- 
ers) of diffusion in cracks [13, Il07l | . Intuitively, it consists 
in the Brownian diffusion of a particle trapped in a ran- 
dom fractal set (a "crack" ) whose pattern resembles the 
graph of a Brownian motion. In the interpretation of 
Eq. (|83p in the crack model, large increments of B\ over 
short intervals correspond to small-width spots along the 
crack, which appear with regularity. 

More generally, time fractional derivatives are involved 
in diffus i on e quations describing transport on fractals [HD, 
0, Il08l. Il09ll . Both diffusion and the medium are then 



irregular. 
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F. Other higher-order diffusion equations 



Without the source term, Eq. (|84[) reduces to the pro- 
totypical higher-order hyperbolic equation 



(d a -vi)u = o. 



(86) 



In contrast to hyperbolic diffusion equations jl05l Il06| , 
parabolic equations such as 



{d a + V A x )u = 



(87) 



describe unconventional diffusion processes^ At first, 
such a differential equation seems unphysical: changing 
the relative sign between spatial and time derivatives is 
tantamount to asking the ink in a bottle of water to con- 
dense back to a single drop. The solution u is not pos- 
itive definite and cannot be interpreted as a probability 
density. These equations, anyway, come from composite 
processes just as their hyperbolic counterparts. In partic- 
ular, they describe diffusion asso ciated with the original 
formulation of the IBM CEqI (but not with IBM as 
later formulated and presented above) or, in the case of 
odd-order kinetic operators V^™ +1 , with th e com position 
of Brownian motion with stable processes |l06| . 



G. Spectral and walk dimension 

Let us consider fractional spacetimes with measure v a . 
The return probability is defined as the spatial average 
of P, i.e., the trace of the heat kernel per unit volume, 

"+00 



1 f^°° 

V((T):=— dg a (x) P(x,x,a) 

*ot J — oo 

Y r+oo r+oo 

= 7T / dg a (x) / dr a ,(k) f k (a) \e(k,x)\' 

J — oo J —oo 

i i r+oo r+oo 

d D x d D kf k (a), 



where V a := j dg a (x) is a divergent total volume prefac- 
tor. Under the rescaling (|55p . we get 



Via) = Ao~— , 
where A is a numerical constant: 

f + °° d D ~x 

J— OO 



(89) 



A = l, 



f_™d D xv a (x) 



0,7 - = 



-oo jD 



d U k 



(27T) 



D J k 



Ml)- 



(90) 

The a dependence of V comes exclusively from (i) the 
rescaling (|55[) and (ii) the volume prefactor in the de- 
nominator. In particular, it depends on the topological 



Parabolic equatio ns of the type H87H were studied in Refs. l86t 
Il044l0^ . Illd - TlT^I and generalized to higher-order operators V™ 
and several other forms. 



dimension D of position and momentum space. It does 
not depend on (i') the measure weight w(k) in momen- 
tum space, nor on (ii') the relative sign between diffusion 
and kinetic operators (hyperbolic or parabolic diffusion 
equation), nor on (iii') the presence of friction or source 
terms, as one can convince oneself by a direct inspection 
(e.g., |86l . ll04l | and references therein). Consequently, all 
these properties will be inherited by the spectral dimen- 
sion. 

The dimensionless coefficient in Eq. (|89[) . formally in- 
determinate, can be regularized so that A = 1. For sim- 
plicity, we specialize to the isotropic case a M = a. The 
i-dependent integrals must be regularized since they di- 
verge at ±oo. Defining two parameters e = e(e) and 
< e <C 1, we write 



+ 00 



f+™d D xv a (x) 



T(a) 



lim 

s->0+ 



lim 

5^0 + 



+oo 

ax 



30 dx x a 



T(a)- 



D 



Jo di 



J eM dx x° 



T(a+1 



e»(e) 



Assuming that Ig l7 > 0, we can take 

The arbitrariness of this procedure does not affect the 
physics, which is encoded in the spectral dimension 



(91) 



da := -2 



leading to the final result 

ds = 



d\nV{a) 
din a 



da . 



(92) 



(93) 



We can distinguish three cases: 

(1) For f) = 7 (including normal diffusion, /3 = 7 = 1), 
ds = dn, whether it be realized by integer or fractional 
differential operators d@ and K.p, a - The important point 
is that the order of the diffusion operator is the natural 
one, i.e., half that of the Laplacian. is a fractal. 

(2) When j3 < 7, diffusion is anomalous and the space- 
time M.® can again be regarded as a fractal. This hap- 
pens, in particular, for integer-order Laplacian and frac- 
tional diffusion (j3 < 1 is assumed). 

(3) As remarked in Ref. (27|, the superdiffusion case 
1 > f3 > 7 does not correspond to a fractal, since 
ds > dn- In particular, for integer diffusion (f3 = 1) 
this is a Levy process. When (3 = 1 and 7 = a, the spec- 
tral dimension coincides with the topological dimension 
of space, ds = D. 

In the case 7 = 0;, the results are summarized in Table 
U for ds and the walk dimension 



'd s 



(94) 
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The walk dimension is also defined through the scaling 



of the second moment, (r ) ~ a 



2/d v 



By a purely di- 



mensional argument, one can infer that 2/dyy = (8/7, 
in agreement with Eqs. {79]) and 1(51)). When d w > 2, 
<iw = 2, d-w < 2, and dw = 1 the process is, respectively, 
subdiffusivc (7 > /?), normal (7 = /3 = 1), superdiffusivc 
(7 < /3), and ballistic (27 = 0). 



TABLE I. Spectral dimension ds and walk dimension dw of 
spacetime for different harmonic structures (Laplacians 
and diffusion equations). M v is fractal only if dw > 2. 





Laplacian 










d s = d H d w = 2 


ds = D > d H d w = a < 2 




d s = /3d H < d H d w = | > 2 


d s = |d H d w = 2f 



A final caveat should be stressed. In the rcgulariza- 
tion procedure, the choice (p?Tj) was possible under the 
assumtpion that T a ,-y > 0. Otherwise, not only the limit 
of the regulator but also the return probability could be 
ill defined, taking negative or complex values. In turn, 
the positive dcfinitcncss of the probability depends on 
the Hcrmiticity of the Laplacian. We have already seen 
that there exists a parameter range such that I Qi7 > 
for the "Riesz choice" of the coefficients in Eq. (fT5"j) . 



V. MULTISCALE PROCESSES AND 
SPACETIMES 

In quantum-gravity theories, the spectral dimension is 
found to depend on the physical length scale £ one is 
probing. In the previous section we found the spectral 
dimension of a fractal model with no scale dependence, 
so that ds is constant. One needs to generalize the dis- 
cussion to a multiscale geometry. The latter will be mul- 
tifractal only in certain ranges of parameter space, where 
ds < <fa- So me basic properties of multiscale measures 
and systems 1 13l4l 18| | were reviewed in Ref. [28j . 

Before dealing with multifractional spacctimes, we 
clarify the problem of multiscaling at the level of the 
diffusion equation, with a = 1 and D = 1. The general- 
ization of Eq. (|54|) to a multiscale process (and a source) 
is simply achieved by summing over all possible values of 
P or 7: 
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u(x — x', a) = S(x — x', a) . 



(95) 

where £^ and C7 ar c dimensionful couplings which de- 
pend on the characteristic scales of the system. The 
parameters /3 and 7 can be also let vary cont inuo usly, 
thus havin g integrations instead of sums [H, Ill9| . or 
both [120]. In particular, by making use of_the so- 
called distributcd-ordcr fractional derivatives 



one can construc t and analyze mult iscale fractional dif- 
fusion equations [HI fl27l fl29Hl33 | . Here we limit the 
discussion to a simple case where the sums in Eq. (|95|) 
are replaced by an integral over a length parameter £, 



d( 



- CiW™] u(x-x\ a, £) = S{x-x', a) , 

(96) 

where the length £ is not identified with (a power of) the 
diffusion parameter a and the range of £ may be chosen 
appropriately (for instance, £ € [£q, +00) in the presence 
of a cut-off). To get the solution of Eq. (|§6"|) for 5 = 0, 
it is sufficient to solve for the integrand, 



d^-a^ ] ] u(x-x',aj) = 0, 



(97) 



where £(£) has been absorbed into From Eq. 193 

the spectral dimension would be (da = 1 here) 



7W 



(98) 



The sum representation is more convenient unless one 
has an argument to choose the profiles (3(£) and j(£) in a 
specific way (as we shall do in Sec. IV C|) . Also, typically 
there is o nly a finite number of terms in physical systems 
HHEl. When the sums are finite, Eq. ([§5"]) docs admit 
analytic solutions in several cases. 

It is important to stress that the number of scales in 
the system determines the number of plateaux or asymp- 
totic regimes in the profile of the spectral dimension ds ■ 
This number, corresponding to the number of Holder ex- 
ponents in a self-similar measure, is given by 

(no. of asymptotic regimes of ds) 
= (no. of /3's) + (no. of 7's) — 1 
= (no. of characteristic scales of the system) + 1 . 

(99) 

Take, for instance, the multifractional-time diffusion 
equation with N diffusion operators and the ordinary 
Laplacian (7 = 1). Then, the number of characteristic 
regimes of the multifractal spectrum is just the number of 
/3's. Assume that a has dimension of a squared length, 
so that the couplings ^ can be written in terms of N 
length scales £ n , n = 1, . . . , N: 



If the ATh operator has ft 
governed by 



£f "~ 2 . (100) 
= 1, the diffusion process is 



N-l 

E 

n=l 



(n 



2(/3„ 



ft 



0, (101) 



2 N is dimensionless. Later we shall better 



where a := a 
specify the scale dependence of the coefficients. 

The discrepancy between the number of fundamental 
scales and the number of regimes is due to the fact that a 
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multiscale phenomenon is always denned by the relative 
size of the scales, i.e., by a hierarchy. This means that 
we can choose any of the N scales £ n to represent the 
variable scale £ probed by a measurement. Thus, there 
are N — 1 (not N) scales with the physical meaning of 
characteristic lengths. The spectral dimension is fixed 
when N = 1; for N = 2 it has two asymptotic values dg ~ 
dsi 2 (with a monotonic0 transient phase in between) in 
the regimes £2 <C l\ and £2 ^ £1', for N — 3 there will 
be two fundamental scales and a plateau ds ~ <is3 at 
intermediate scales; and so on. 

We need only one scale to obtain a multifractal or, 
more generally, a geometry with scale-dependent spec- 
tral dimension. Therefore, as a first approximation, it is 
sufficient to study diffusion equations with two diffusion 
operators d^ 1 and d@ 2 or with two Laplacians d 2 ^ 1 and 

®\x\ ■ Certain realizations of QEG can be regarded as a 
two-scale system, where there is an intermediate genuine 
plateau [361 ] . 



A. Multiscale Levy process 

The first example is an interaction of Gaussian and 
anomalous dynamic s which can describe certain turbu- 
lent media [HH 1 1251 ] . The diffusion equation is 

(d„ -d 2 x - C^g) u{x -x f ,a) = 0, (102a) 



u(x — x , 0) = S(x — x) . 



< 7 < 1 , (102b) 



where [a] — —2 and £ is a coupling constant which, by 
a dimensional argument, can be written in terms of a 
characteristic length scale £*: 

C = ^ 2(1 ~ 7) , K1=2(l- 7 ). (103) 

The Fourier transform of the spatial part is 

(fc 2 + C|fc| 27 )u(fc,a) ^ £- 2 [(^k) 2 + \£^}u(k 7 a) , (104) 

stating that the transport is normal at small scales fc _1 <C 
£* and of Levy type at large scales hr 1 ^> I*. At scales 
k^ 1 ~ I* the two behaviors compete. The normalized 
analytic solution can be written in different ways [68j . 
For us, the following is convenient: 



+00 



u(x-x',a) = ^2(-l) n y n (a)\x-x' 



1 



<2n)\TT 



+00 



dk k 



(105a) 
(105b) 



9 A local extremum, for instance a minimum value dg min at some 
point < f m in < +00 would signal the presence of another scale 
because this feature could not be removed by a finite conformal 
rescaling of diffusion time. 
10 Some remarks on the case with a finite but arbitrary number N 
of diffusion operators are given in Ref. [5011 . 



where all barred quantities are dimensionless: 

k : ^ £-± k . 



x 



(106) 



The coefficients y n can be expanded as a series for large 
and small £^ 2 cr. In particular, we are interested in the 
return probability and for x — x' only the n = term 
contributes: 



V{a) = y (a) 



(107) 



Therefore, the spectral dimension is ds ~ 1 at small 
scales and dg ~ I/7 > 1 at large scales. This is con- 
sistent with Eq. (|98j) with (3(£) = 1 and a profile 



7W 



1 , £ < 4 
7 < 1 , £ > U 



(108) 



From the perspective of quantum spacetimes, this 
model is multiscale but not multifractal, since ds > dn. 
Profiles of ds over-shooting the Hausdorff and topologi- 
cal dimensions of space appear also in lattice-based [1 51 ] 
and noncommutative geometries (22| . Due to the pres- 
ence of a characteristic scale (the noncommutative fun- 
damental scale [13], or the lattice cell siz4]3 [HJ, or the 
label-dependent length of the edges of a labeled graph 
[l5| ) , the geometric information in the diffusion equation 
determines a nontrivial spatial generator. On the other 
hand, the diffusion operator is assumed to be the integer 
one d a , so these models roughly mimic certain properties 
of Levy processes. In these cases, the flow is not mono- 
tonic since the scale acts as a minimum length cutoff, so 
ds —> in the UV. Since ds — > D at large scales, there 
is one local maximum in ds(£) (for labeled graphs there 
are, in fact, several extrema in ds fl5|). 



B. Fractional telegraph process 

A fractional diffusion equation with multiple diffusion 
operators d@ n admits a neat stochastic interpretation in 
the case N = 2. For the purpose, w e recall s ome results 
on the so-called telegraph processes |l35l4l38| . After this 
short review, we shall make the connection with multi- 
scale quantum spacetimes apparent. 

A telegraph process V in time a (here [a] = —1) is 
defined as 



V{a) = V{0) (-l)^( ff ) 



(109) 



where V(a) is the velocity of a particle at time a run- 
ning on the real line, V(0) is the initial velocity which is 



11 In the case where the underlying lattice- type construction is 
taken as a rcgularization, these are lattice artifacts. 
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±c with equal probability, and N is the cumulative num- 
ber of events of a homogeneous Poisson process. Thus, 
the velocity of the particle flips direction at times obey- 
ing a Poisson distribution, hence the name "telegraph." 
The position of the particle at time a is the integrated 
telegraph process 



T(a)=V(0) / ds(-l)^ 



(110) 



and i ts pr obability distribution obeys the diffusion equa- 
tion [he mi 



(dl + 2Xd a -c 2 dl)u(x-x\cr) = 0, (111a) 
u(x - x', 0) = 5(x - x') , (111b) 

where A > is the rate of the Poisson process. The vari- 
ance of the process can be shown to scale as (r 2 ) ~ a 
at large a. When A, c — > +oo and the ratio A/c 2 re- 
mains constant, Eq. (|111[) reduces to the ordinary dif- 
fusion equation |136(. The limit A — > +oo means that 
changes in the speed, abrupt in the telegraph process, 
take place continuously. 

The composition of a Brownian motion B with an in- 
tegrated telegraph process, 



X 



DBM 



B[T(a)] , 
iB[-T(a)] , 



T(a)>0, 
T(cr) < 0, 



is called delayed Brownian motion (DBM) |l06| . and its 
probability distribution is governed by the quartic diffu- 
sion equation 



dl + 2Xd a - —dt ) u(x - x', a) 



u(x — .t', 0) = S(x — x') . 



A/c 2 



, (112a) 
(112b) 

const, the de- 



In the double limit A, c — > H 
laycd Brownian motion reduces to Eq. (|87|) . Other combi- 
nations of telegraph and Brownian processes ar e possible , 
leading to different quartic diffusion equations [l05l Il06j . 

Finally, one can devise the fractional analog of the tele- 
graph equation [85l . l98j : 

(df + 2A^ - c 2 8l) u 20 (x - x',a) = , (113) 



with initial condition given by Eq. (|67[) or (JHHJ) . When 
A = 0, Eq. (|113|) reduces to the fractional diffusion equa- 
tion (|66[) . while for (3 = 1 it is the telegraph equa- 
tion . Ge neralizations of Eq. (|113p were inspected in Refs. 

mini. 

The analytic solutions of the telegraph equation (|lllj) . 
of the delayed-Brownian- motion equation (|112[) . and of 
the fractional teleg raph equ a tion (|113p were found, re- 
spectively, in Refs. [l3l[l38j], [l06|, and [H|. The latter, 
which is non-negative and unique, is given as a Fourier 



integral: 

2 r+oo 

u 2f }(x-x',a) = — / dku 2 p{k,o)e- lkx , (114) 



W2/3(fc,cr) = - 



-oo 

V- 



A + T]- 



X + V+ 

where E is the Mittag-Lefllcr function (|22l) and 



V± 



±yj\ 2 - c 2 k 2 - A. 



(115) 



(116) 



The variance of the process scales as (r 2 ) ~ cr' 3 in the 
limit of large a. 

In the special case (3 = 1/2, the fractional telegraph 
equation becomes 



(d a + 2X83 - c 2 9 2 ) m(x -x',a) = 0. 



(117) 



Its solution can be written explicitly as a one-parameter 
integral: 



ui(x — x\ a) 



1 



-As 



ds ■ 



S(r - s) + S(r + s) 



+8(s-r) XI (a^s 2 - r 2 ) 
d s I (Ws 2 -r 2 )]} , 



(118) 

where r = \x - x'\/c and I (z) = E^ToO / ' 2 f n I '(nl) 2 is 
the modified Besscl function of the first kind and zcroth 
order. 

In the double limit A, c — >• +oo, A/c 2 — s- const, the 
process underlying (|117p reduces to an iterated Brownian 
motion and, in fact, the limit of Eq. (|118p is ([55]). The 
probability distribution u\ is associated with a telegraph 
process with Brownian time, 



JW(<r) =T[\B{a) 



(119) 



which we call also fractional telegraph process. The com- 
posite process Xftp describes the random motion of a 
particle during a time interval of length \B\, so that at 
time a the particle is located in the random spatial in- 
terval (~a\B(a)\, a\B(a)\) (98|. This motion is governed 
by the multifractional diffusion equation (|117l) . Due also 
to the different orders of composition of T and B, the 
fractional telegraph process Xftp is different from the 
delayed Brownian motion Xdbm- 

We now end the digression into probability theory and 
recast these results in the language of multiscale space- 
times, as follows. We have to reshuffle the units of the 
variables and constants to agree with our past notation, 
where [cr] = —2. Setting c = 1 and I ~ 1/(2A) as the 
probed scale, Eq. (|117|) becomes 



d s + 4 9f - ildl ) Ul = , 



(120) 
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FIG. 4. The spectral dimension ds(i) of a D = 1 multifractal 
space whose stochastic properties are given by a fractional 
telegraph process. 



where, as before, a — £~ 2 a is dimensionless. In the limit 
I 3> £*, diffusion in spacetime is Gaussian and described 
by a Brownian process. At small scales f « 4, on the 
other hand, one reaches a regime where diffusion is frac- 
tional and described by an iterated Brownian motion. In 
between, diffusion in quantum spacetime obeys the law 
of a fractional telegraph process. 
The return probability 



Via) =u 1 (0,a) 



1 



1 



1 

4? 



+ 00 



Jo 



2( 



Ji 



ds t 

s 

Yi 



(121) 



can be manipulated to give a sum of generalized hyperge- 
ometric functions with argument a j£ 2 — (£*/£) 2 a. From 
that, we can obtain the profile of the spectral dimension 
ds(£)t which is shown in Fig. [4] (the portion of the curve 
below ds = 0.65 down to ds = 1/2 is not shown because 
of bad convergence of the code). As expected, 



ds 



(IR), 
(UV). 



(122) 



The results above can be generalized to D dimensions 
because the Fourier integral (|114p can be split into an 
angular and a radial integral, d D k = dVln-idk fc D_1 . 



C. Multifractional measure 

All the above examples are realizations of multifrac- 
tional spacetimes where a = 1 and the Hausdorff dimen- 
sion of position space equals the topological dimension, 
dji = D. The spectral dimension is anomalous cither 
because of a fractional diffusion operator or due to a 
fractional Laplacian, or because of both. Inclusion of 



a nontrivial measure weight for position space makes g?h 
anomalous, too. Formally, the extension of fractional 
to multifractional spacetimes is simply achieved by sum- 
ming (or integrating) over all possible values of a (26l.[28|. 
A multifractal action with a finite number of charges a 
thus reads 



S 



N 
n=l 



dP X V an (x) £ Q „, 



(123) 



where g n are some (dimensionful) couplings and £ Qi7 is 
the Lagrangian density, possibly dependent on a and/or 
the order of the kinetic operator. The aim is to find the 
scale-dependent spectral dimension ds- This entails a 
modification of the diffusion equation and the introduc- 
tion of scale dependence in the parameters a, (3, and 7. 
Due to its phcnomcnological character, we do not expect 
to obtain a unique multiscale extension. However, one 
can restrict the possibilities by an educated guess. 

The general diffusion equation on the spaces described 
by (fT23|) is 



N 

E 



{tnd^ - C»£® ) P(x,x',a) = S{x,x',a) , (124) 



given some initial condition P(x, x', 0). If we fix j n = 7 
and order the N — 1 scales of the system as £\ < £2 < 
■ ■ ■ < £n-1i we can argue that the coefficients £ ra have 
the forms = 1 and 



(125a) 



( 



2l 



ti = 2,...,JV- 1 



(125b) 



where I = £n > £n-i- To show this, we notice that the 
Laplacians all have the same order 27, so the coefficients 
Q n all have the same scaling dimension. By fixing the 
scaling of the suitably, we can always make Q n dimen- 
sionless. This means, in particular, that we can write 
as the ratio of some length scales, ( n = (lA,n/lB,n) 9 ■ 
Without loss of generality, one can choose q = 27 so that 
the spatial generator of the diffusion equation can be ren- 
dered dimensionless, in the form J2 n (^A,n) 2 ' 1 l^-f,a n - Now, 
the nth term dominates over the others at scales I <C £ n , 
so we could set Ia,u = £n and, tentatively, Is.n = £■ How- 
ever, at scales smaller than £ n -i the (n— l)th term takes 
the lead, so the smallest possible scale £ at which the nth 
term dominates is £ ~ l n -\- Therefore, the correct choice 
is Ib,u — £ — £n-i- In other words, the dimensional flow 
is always measured starting from the lowest of two scales 
l n -\ to the next £ n , and relatively to the latter, which 
sets a gauge for the rods. Beyond the smallest scale £\ 
there is nothing else to compare with and £q — 0. Since 
£ = £n is the probed scales, (n = 1 by definition. 
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In the Gaussian case /3 = 7 = 1 and just two entries 
(N = 2, c*2 = 1, oti = a*), dimensional flow is such that 



£> , £ » 4 (IR) , 

L>a* = 2 , I < 4 (UV) , 



(126) 



with no intermediate regimes in between. Explicitly, Eq. 
(|124j) becomes 



P(x,x ,a) = S(x, x , a) . 



(127) 

where we set (2 = 1 and £i = £ = (£*/£) 2 according to 
Eq. (|T25|) . When 5 = 0, the UV and IR asymptotics of 
the solution arc obvious. 

Equations (|124[) and ()127[) are not much prone to ma- 
nipulation, since fractional momentum transforms are 
not easi ly g eneralizable to the multifractional case when 
D > 2 (29J. However, we can reinterpret the sum over 
a by imagining the dimensionless parameter a = a(£) 
to depend on the probed scale. Equation (|123p can be 
recast as 

r+00 /• 
S = j^ d£g(£)J d D xv aW (x)C aWMe) , (128) 

while the diffusion equation (|124[) generalizes Eq. (|M|) to 

J d£ [t(£)9£ W - CWmmq] P(x,x',a,£) 

= S(x-x',a). (129) 

In the absence of source, it is sufficient to solve for the 
integrand, 



d m _ {(£)K* {e)M J P(x, x',a, 0=0, (130) 



where, as before, we absorbed £(£) in ((£). All the cal- 
culations of the case with fixed dimensionality are trans- 
posed with the replacement (a,/3,7) — > (ct{£), P{£ ), j(£ )) 
and the spectral dimension is the generalization of Eq. 



ds(£) 



i{£) 



da(£), d H (l) = Da(e). 



(131) 



We can plot the spectral dimension for any given pro- 
files a(£),j3(£),^(£). To obtain a multifractal, ,9/7 < 1 
throughout the whole evolution, so for simplicity we as- 
sume 7 = P = 1. In Ref. [28j], several Ansdtze were given 
for an a(£) with the desired asymptotic behavior (|126[) F^1 
They are all monotonic and lead to the same plot qual- 
itatively, the only change being in the time-scale of the 
process. Here, however, we choose a profile which can be 



motivated as a realistic approximation of the sum (| 1241) 
with 7„ = 1 for all n. Consider first the N = 2 case with 
Q'l = a* and Qt2 = 1. In one dimension, 



(1 + C1) 



1 



Ci (1 - a*)(3 - a*) 



1 + Ci«« 
1 + C1 

a* 



-9s 

X 



where 



ax{£) 



1 + Ci 4a; 2 

(i + Ci)/cf iQl(£) + 



1 + CiW a* 



P 



Ci (1-a*) 2 



1 + Ci 



4x 2 



P, (132) 



1 + C1W 



(133) 



For both small and large £1 the kinetic term in Eq. 
(|132j) dominates over the potential term, so the profile 
(|135p defines an effective fractional charge a c g ai(^) 
throughout the dimensional flow. With N coefficients 
a ni Oi 7y = 1, the effective fractional charge reads 



a N -\{£) ■-- 



(134) 

In fact, this is nothing but the average (a) of the coeffi- 
cients a n with respect to the weights £„. 

The spectral dimension for D = 4 in the one-scale case 
(N = 2, a* = 1/2) is shown in Fig. [5ja), in agreement 
with|U A two-scale profile with ot\ = 1/2, «2 = 1/3 and 
h = 104 is plotted in Fig.ETb): 



M*) := 



I + K1W + K2W 



c - if 



I + C1W + C2W 
2 

c 2 = 



(135a) 
(135b) 



At £ = (beginning of the flow, UV critical point), da, = 2 
in four dimensions. At £ ~ £1, the spectral dimension 
acquires the minimum value dg = 4/3. At scales £ < £ 2 , 
the d iffusi on process corresponds to a recurrent random 
walk jl40| . where dw > <fe (ds < 2) and each site of the 
walk within a given radius is visited several times. Well 
above the larger critical scale, £ ^ £2, both dimensions 
hit the IR value ~ 4: 



D, 
Da 2 = %, 

D 
2 ' 



Dai 



£^>£ 2 ->£i (IR) , 
?i ~ £ <C £ 2 (intermediate) , 
£ < £1 < £ 2 (UV) . 

(136) 



12 Equation (1.1) of Ref. [2Sll has a typo and should read ot(i) 
1 + (q» -l)/[l + (^/£) Q *- 1 ]. 



The two-scale spectral dimension of Fig. [5T b) repro- 
duces the dimensional flow of QEG [IJ l30t • 



20 




0.01 0.1 1 10 100 ' 0.01 0.1 1 10 100 1000 

tit. f/t, 



FIG. 5. The spectral dimension ds(l) in D — 4 for a multifractional model and normal diffusion (/3 = 1 = 7, ds = du) in four 
dimensions, with profile (|133[) (a) and (|135p (b). 



VI. DISCUSSION 

Dimensional flow, the change of spacctime dimen- 
sionality with the scale, is a trademark of quantum- 
gravity scenarios which appears in different guises but 
bears a simple, limited set of general characteristics. 
We have listed and classified these characteristics in a 
portable fashion (i.e., independently of the quantum- 
gravity model) according to their degree of dependency 
on the details of the diffusion equation. While the num- 
ber of asymptotic regimes (plateaux) for the spectral di- 
mension, their values and positions are physical quanti- 
ties related to the hierarchy of fundamental scales of the 
system, the intermediate, transient regimes connecting 
them depend on the specific realization of the diffusion 
equation. We have dissected the latter in its constituent 
elements and classified a number of ways to introduce a 
multiscale structure in it. While doing so, we pointed 
out how uncharted territory in quantum gravity can be 
explored using the maps of other branches of physics 
and mathematics such as diffusion and stochastic the- 
ory. Anomalous and multiscale diffusion equations, in 
fact, have been known since long ago in these disciplines. 
Some of the contributions of the present paper amount to 
link composite stochastic processes to multiscale space- 
time geometries, as in the case of the fractional telegraph 
process (even in probability-theory literature, to the best 
of our knowledge, the latter seems not to be directly asso- 
ciated with a multiscale system) , to "locally" characterize 
asymptotic regimes in the dimensional flow with specific 
stochastic processes, and to provide a qualitative (Sec. 
HI)) and quantitative (Sec. [Vj general analysis of dimen- 
sional flow. In parallel, we have improved the status of 
knowledge of multifractional spacetimes regarding frac- 
tional Laplacians, diffusion, and the analytic treatment 
of multiscale configurations. 

The single-scale example is clear-cut in the way it il- 
lustrates how to associate a given dimensional flow or 
portions of it with specific types of stochastic processes. 
Suppose we have a quantum-gravity model with space- 
time spectral dimension following the monotonic profile 



of Fig. Q] Different diffusion equations can give the same 
profile qualitatively, with the same asymptotics but a dif- 
ferent slope in between [Figs. 2] and 03a) and, e.g., Fig. 1 
of Ref. [4l|] . As we have seen, one such diffusion equation 
governs what is known as a fractional telegraph process, 
which is a telegraph process with Brownian time (Sec. 
IV Bl) . In this stochastic process, the test particle expe- 
riences abrupt changes of speed at times governed by a 
Poisson law, in turn parametrized by a Wiener process 
(i.e., Brownian motion). At large scales, towards the up- 
per right plateau of the figure, this process reduces to 
an ordinary Brownian motion. At small scales, in the 
lower left plateau, the process reduces to another com- 
posite process, called iterated Brownian motion, which is 
a Brownian motion with Brownian time. Iterated Brown- 
ian motion describes the diffusion of a Brownian particle 
in a fractal medium. Now, suppose the quantum-gravity 
model at hand does not realize this particular diffusion 
equation (for instance, multifractional theory). Since di- 
mensional flow is essentially the same, we can still in- 
terpret the one-scale anomalous diffusion equation as de- 
scribing the diffusion of a particle in an irregular fractal 
medium, although the process is not really a fractional 
telegraph process. Nevertheless, the existence of a spe- 
cific stochastic process in the same "equivalence class" of 
flows (in the example, the class of monotonic single-scale 
flows) gives a sharper characterization of the physics un- 
derlying diffusion in an anomalous spacetime. This is 
also true asymptotically; in the example, the asymptotic 
class of flows leading to ds ~ D/2 is that of the IBM. 



Through the exact and approximate solution of the 
various diffusion equations, we have replaced spotwise 
knowledge of asymptotic values of the spectral dimension 
with a continuous analytic control over dimensional flow 
in a model-independent fashion. In particular, control 
over multifractional spacetimes has been ameliorated. 
This should open up further advancements in other as- 
pects of the field as well as its application as an effective 
framework to other theories of quantum gravity. 
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